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ABSTRACT 

For  an  effective  adhesively  bonded  reinforcement  of  a  (typically  cracked)  plate, 
sufficient  load  must  be  transferred  by  the  adhesive  into  the  reinforcement  to  prevent 
the  underlying  damage  from  growing.  Under  severe  load,  the  adhesive  may  yield 
plastically.  Notwithstanding  possible  adhesive  failure,  this  may  be  beneficial  to  the 
reinforcement  by  reducing  the  peak  stress  adjacent  to  the  crack.  In  this  paper, 
characterisation  of  this  stress  reduction  compared  to  the  (non-yielding)  elastic  case, 
was  sought  by  examining  the  influences  of  configurational  parameters  including  plate, 
adhesive  and  reinforcement  moduli,  and  adhesive  yield  stress.  Finite  element  (FE) 
analyses  were  conducted  for  a  two-dimensional  section  through  a  double-sided 
(s}rmmetric)  lap  joint,  representative  of  a  typical  repair.  Stress  reductions  in  the 
reinforcement  of  the  order  of  25%  were  foimd.  The  adhesive  yield  was  shown  to  be 
dominated  by  shear  stress,  and  thus  the  adhesive  behaved  essentially  one- 
dimensionaUy.  The  linear  increase  in  plastic  zone  length  with  applied  load,  as 
predicted  by  the  Hart-Smith  one-dimensional  theory,  was  in  good  agreement  with  the 
FE  results.  However  the  observed  load  transfer  length  was  6-18%  longer  than 
predicted. 
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Finite  Element  Analysis  of  the  Double  Lap  Joint 
with  an  Elastic-Plastic  Adhesive 


Executive  Summary 

At  least  three  methods  are  available  to  restore  serviceability  of  an  aircraft  with  damage 
such  as  cracking  in  structural  components  and  battle  damage.  These  include  replacing 
the  damaged  part,  riveting  or  bolting  on  a  reinforcement  plate,  and  application  of  an 
adhesively  bonded  reinforcing  patch  over  the  damage.  This  third  medrod  features 
reduced  downtime,  and  avoids  further  damage  caused  by  driUing  more  holes.  It 
requires  that  the  reinforcement  be  stiff  enough  to  take  sufficient  load  to  prevent  the 
underlying  damage  (crack)  from  growing,  but  not  so  stiff  that  it  substantially  alters  the 
stress  distribution  originally  designed  for  in  the  structure.  This  work  was  undertaken 
to  address  two  issues  concerning  bonded  repair,  the  first  being  the  accuracy  of  a  RAAF 
engineering  standard  covering  this  type  of  repair.  The  second  issue  is  the  high  level  of 
stress  in  the  reinforcement,  and  whether  this  can  be  reduced  enough  to  eliminate 
possible  failure  of  the  reinforcement. 

The  RAAF  engineering  standard  is  based  on  a  simple  model.  This  paper  seeks  to 
establish  a  basis  for  the  accuracy  of  the  model,  or  to  produce  more  accurate  formulae. 
The  focus  is  on  a  (symmetric)  double-sided  lap  joint,  representative  of  the  worst-case 
repair  of  a  through-cracked  plate  but  without  bending.  There  is  a  brief  examination  of 
one-sided  repairs  where  peel  stresses  and  bending  are  more  important.  Load  transfer 
from  plate  to  reinforcement  occurs  via  (principally  shear)  stress  in  the  adhesive.  Finite 
element  (FE)  simulations  using  the  PAFEC  computer  package  produce  the  underlying 
data  for  analysis.  With  tiie  symmetric  lap  joint,  the  dominant  adhesive  stress  is  found 
to  be  shear.  This  explains  the  accuracy  of  tiie  one-dimensional  model  which  only 
allows  for  longitudinal  stresses  in  the  plate  and  reinforcement,  and  shear  in  the 
adhesive. 

This  work  extends  the  model  to  include  the  situation  where  the  adhesive  yields 
plastically.  Above  a  specified  (yield)  stress,  the  adhesive  distorts  further  without 
canying  any  higher  stress.  Although  this  yielding  is  normally  dependent  on  all  the 
stress  components,  the  dominance  of  shear  meant  the  adhesive  was  again  well 
described  by  a  1-dimensional  model.  This  yielding  extends  the  lengtii  over  which  the 
adhesive  transfers  load  from  the  plate  to  &e  reinforcement.  If  the  adhesive  does  not 
fail,  this  has  the  beneficial  effect  of  lowering  the  peak  stress  occurring  in  the 
reinforcement  adjacent  to  the  crack.  This  reduction  is  significant:  from  around  1.6 
times  the  average  stress,  to  1.3  times.  Characterisation  of  this  beneficial  effect  is  sought 
by  examining  its  dependence  on  configurational  parameters  including  plate,  adhesive 
and  reinforcement  moduli,  and  adhesive  yield  stress.  The  linear  increase  in  plastic 
zone  lengtih  with  applied  load,  as  predicted  by  the  model,  is  closely  simulated  by  the 
FE  analysis  although  the  observed  load  transfer  length  is  6-18%  longer  than  calculated. 

This  work  has  shown  that  the  simple  model  imderlymg  the  RAAF  design  standard  on 
bonded  repair  accurately  describes  the  behaviour  of  the  repair  components  under  load. 
The  model  was  extended  to  allow  for  plastic  yielding  in  the  adhesive.  This  was  shown 
to  have  the  benefit  of  lowering  the  peak  stress  occurring  in  the  reinforcement. 
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1.  Introduction 

In  the  contexts  of  stractural  integrity  and  life  extension,  adhesively  bonded 
(reinforcement)  patches  have  been  used  to  repair  cracked  metallic  plates  (Baker  and 
Jones(l),  Rose  et.  al{2)).  For  an  effective  repair,  the  adhesive  must  transfer  sufficient 
load  into  the  reinforcement,  to  prevent  the  underlying  damage  (crack)  from  growing. 
At  the  same  time,  the  reinforcement  shordd  not  be  so  stiff  that  it  substantially  alters  the 
stress  distribution  originally  designed  for,  in  the  structure. 

In  severe  cases,  the  adhesive  may  yield  plastically.  If  the  adhesive  does  not  actually 
fad,  this  plastic  yielding  may  benefit  the  repair  by  lowering  the  peak  stresses  occurring 
in  the  reinforcement  adjacent  to  the  crack.  However,  it  is  desirable  that  the  yielded 
zone  be  small  enough  to  preclude  subsequent  adhesive  degradation. 

This  work  focusses  on  the  double  lap  (symmetric)  joint,  as  a  model  for  two-sided 
repairs  as  indicated  below.  The  aim  is  to  seek  design  parameters  that  characterise  the 
repair  effectiveness,  with  greater  accuracy  than  that  obtained  with  the  simple  Hart- 
Smith(3)  one-dimensional  model. 

Peel  stresses  and  bending  are  known  to  be  important  in  adhesive  fadure,  partiadarly 
for  a  one-sided  repair.  These  are  briefly  examined  but  are  not  the  focus  of  this  analysis. 
Reinforcements  are  usuady  tapered  towards  the  edges  to  minimise  these  effects,  but 
this  may  be  restricted  in  situations  where  the  patch  size  is  limited. 

In  this  work  finite  element  analyses  using  the  PAFEC  level  8.1  package  were 
conducted.  These  results  were  initiady  benchmarked  against  known  analytical  residts 
for  an  elastic  adhesive.  They  were  then  used  to  produce  data  for  the  case  where  the 
adhesive  becomes  plastic,  yielding  above  a  nominated  von  Mises  (Henshed(4))  yield 
stress. 


1.1  The  Geometry  of  the  Simplified  Problem 

The  bonded-reinforcement  repair  of  a  crack  is  a  three-dimensional  problem,  as  shown 
in  Figure  1.  In  the  present  analysis,  a  two-dimensional  sHce  (A-B  in  Figure  1)  is  taken 
through  the  plate  and  reinforcements  including  the  critical  region  in  the  vicinity  of  the 
crack.  This  models  the  repair  as  a  double  (symmetric)  lap  joint,  representing  the  "worst 
case"  repair  where  the  crack  extends  right  across  the  plate.  In  view  of  the  symmetry 
either  side  of  A-B,  there  is  no  displacement  out  of  the  A-B  plane.  This  constrains  the 
problem  in  Figure  lb  to  plane  strain. 
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(a) 


(b) 


Figure  1.  (a)  A  reinforcement  patch  repair  and  (b)  the  corresponding  double  lap  joint  with 
(symmetry)  reduced  region  shown  dashed. 


Two  further  assumptions  permit  the  reduction  to  a  conceptually  and  analytically 
simpler  problem  shown  by  the  dashed  rectangle  in  Figure  lb. 

1.  The  edge  of  the  reinforcement  and  edge  taper  are  far  away  compared  to  the  load 
transfer  length  across  the  adhesive.  This  leads  to  a  region  far  from  tihe  crack,  and 
edge  taper,  of  constant  and  equal  strain  Sxx  in  the  plate,  adhesive  and  reinforcement. 
The  problem  can  then  be  tnmcated,  replacing  the  part  to  the  right  with  uniform 
stresses  across  the  plate,  adhesive  and  reinforcement  in  proportion  to  tiieir  moduli. 

2.  Symmetry  restricts  the  problem,  to  a  region  above  the  centre  tine  and  to  the  right  of 
the  crack  with  appropriate  constraints  as  shown  in  Figure  2. 


Mjc(0,y)=o 

4 

,  Reinforcement 

♦ 

^C5cc=^/1000 

^  C5cc  =^/1000 

\ _ at _ 3 

1 

4 

Adhesive  ^ 

Plate 

c^=£^/1000 

0  ^  Uy  (x:,0)=0  (none) 

(m3;(0.02,0)=0) 


Figure!.  Reduced  double  lap  joint  problem  showing  the  configurational  parameters  of 
thickness  (t),  Young's  modulus  (E)  and  constraints:  displacement  (u)  and  the 
applied  stresses  for  the  reference  load  (constraints  for  the  one-sided  repair  in 
parentheses). 
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For  this  problem,  the  centre  line  (i/=0)  constraint  is  that  the  vertical  displacement  Uy=0. 
The  constraints  along  x=0  are  that  the  crack  face  (plate)  has  oi*=cr*y=0  and  the 
reinforcement  Ux=0  along  ;t:=0.  For  the  one-sided  case,  the  Uy=0  constraint  is  eliminated. 
An  additional  constraint  is  then  needed  to  fix  the  overall  vertical  displacement:  say 
Wy(0,0)=0.  To  simplify  specification  of  the  Ux{0,y)=0  boundary  condition  across  the 
reinforcement,  the  plate  has  been  terminated  at  3:=0.02mm  in  the  unloaded  state.  This 
means  the  crack  has  a  finite  but  small  width  at  zero  load. 

The  simplified  problem  with  boundary  conditions  is  illustrated  in  Figinre  2.  It  is 
assigned  a  depth  into  the  page  of  a;=100mm.  The  stresses  on  the  right  end  of  the  plate 
at  the  100%  reference  load  level  (denoted  by  Fr)  then  convert  to  forces  given  for  the 
plate  by  (Table  1):  F=cr.tP.w=22MN.  The  parameters  chosen  for  the  analysis  are  set  out 
in  Table  1.  Note  that  t?  is  the  half-thickness  of  the  plate  shown  in  Figure  1. 


Table  1.  The  parameters  used  in  the  analysis. 


Parameter 

Plate 

Adhesive 

Reinforcement 

Material 

A1 

FM73 

Boron  fibre  (10  ply) 

2024 

composite 

Young's  Modulus:  E  (GPa) 

72.4 

1.890 

210 

Poisson's  Ratio:  v 

0.33 

0.35 

0.30 

Thickness:  t  (mm) 

3.10 

0.20 

1.27 

Depth:  w  (mm) 

100 

100 

100 

Relative  stiffness:  Et  (mm  MPa) 

224 

0.4 

263 

Stress  apphed  at  100%  load:  cr  (MPa) 

72.4 

1.890 

210 

Force  at  100%  (Fr):  F  (kN) 

22.4 

0.04 

26.3 

Shear  modulus:  //  (GPa)* 

27.2 

0.70 

79.6 

von  Mises  yield  criterion: 

low:  37.0 

<7^  (MPa) 

high:  64.1 

*  The  shear  modulus  is  given  by  //  = 


E 

2(1  + v) 


1.2  Specific  Cases  examined 

The  following  cases  were  examined: 

1.  An  elastic  calculation  for  the  two-sided  (symmetric)  joint  is  the  benchmark  case. 

2.  An  elastic  calculation  for  the  one-sided  reduced  lap  joint  problem  extracted  from 
cracked-plate  configurations  as  indicated  above,  is  compared  to  1. 

3.  The  two-  sided  joint  is  then  examined  for  two  different  plastic  yield  stress  cases. 
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4.  The  one-sided  plastic  case  is  presented  with  its  higher  peel  stresses,  important  for 
failure  in  some  real  repairs. 

5.  Two-sided  cases  with  varying  plate,  adhesive  and  reinforcement  moduli  are 
examined  to  determine  the  effects  of  plasticity  in  these  differing  circumstances. 


2,  Mathematical  theory 

2.1  Separation  of  von  Mises  Stress  Components 

There  are  a  number  of  analytical  formvdae  needed  in  this  analysis.  In  particular,  the 
von  Mises  peld  condition  (Henshell(4))  is  based  on  the  effective  stress  ae  where: 

20-'  =  +  (^yy  ■  (1) 

This  has  been  split  into  the  major  (one  dimensional)  component  <Jxy  and  an  additional 
component  oi. 

^  +  ^xy}-  (2) 

Restricting  the  problem  to  plane  strain,  the  following  equatioirs  are  obtained. 

=  CTyz  =  0.  =  v(<T,,  +  <y^) 

o-'  =  j{[l- v+v'](cr^-cr^^)'+[l-2v]V,,f7^^|  ^  ^ 

Substituting  in  the  value  for  FM73  adhesive,  v<=0.35,  this  reduces  to 

^  "  6  “  ^yyf  •  (4) 


2.2  Load  Transfer  Length  and  Yield  Load:  One  Dimensional  Model 

The  Hart-Smith(3)  one-dimensional  model  allows  for  longitudinal  stresses  only  in  the 
plate  and  reinforcement,  and  shear  in  the  adhesive.  By  considering  equihbrimn  of 
these  forces  acting  on  the  components  (Figure  3a),  equations  may  be  obtained  for  the 
stresses  shown.  The  adhesive  stress,  cr^{x)  is  found  to  decay  exponentially  with  the 
form  exp{-/3x).  The  load  transfer  length,  is  given  in  this  model  by 


E 


Here,  R  and  P  refer  to  the  reinforcement  and  plate  respectively. 


/^\l  +  5 


S  E^t 

—  where  the  stiffness  ratio  is  5  = - ^ 


EU, 


(5) 
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This  Hart-Smith  model  predicts  a  linear  increase  in  plastic  zone  length  with  applied 
load.  This  can  be  seen  in  the  following  figure  showing  the  oiy  stress  in  the  adhesive 
(the  only  stress  considered  in  the  adhesive  in  this  model)  after  some  plastic  yield  with 
the  plastic  zone  length  Ip.  The  total  force  transmitted  across  the  adhesive  (Figure  3b)  is 
given  by  the  area  xmder  the  curve, 


Ipxaly+rerly  exp(-/3[x  -  Ip])  6x 


=  w  X  E  tp  X - - 

1000 


(6) 


The  1000  arises  because  the  reference  load  (Fr)  is  defined  as  that  loading  where  the 
applied  stresses  are  1/1000  of  the  Yoxmg's  modulus  for  the  plate,  adhesive  and 
reinforcement  respectively  (Figure  2).  As  the  applied  load  (F)  increases,  the  plastic 
zone  (Ip)  must  extend  enough  (Ap)  to  transmit  the  additional  load.  This  will  be 
proportional  to  Apxa^xy.  For  a  one  dimensional  analysis,  the  von  Mises  5deld  criterion 
(equation  (3)  withcrfl=0),  gives  a^xy=o^f'i3. 


<7^(x  +  Sx) 
a^„(x  +  Sx) 


Figure  3.  (a)  Infinitesimal  element  used  in  deriving  the  one  dimensional  theory,  showing  the 
stresses  and  displacements  considered,  (b)  Plastic  zone  of  length  Ip  with  exponential 
tail  in  the  adhesive  shear  stress. 


It  is  straight  forward  to  carry  out  the  intergation  in  equation  (6)  and  rearrange  to  show 
the  predicted  linear  increase  in  plastic  zone  lengdi  with  applied  load: 


E'^tp  F/F^  1 

<  ""  1000  p 


p 


\Fy  J 


(7) 


The  second  equality  above  introduces  the  load  Fy  needed  to  initiate  plastic  )delding. 
The  first  shows  that  the  plastic  zone  growth  rate  is  only  dependent  on  Ae  plate 
properties  and  the  adhesive  yield  stress. 
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3.  Numerical  Analysis 

3.1  The  finite  element  mesh 

The  mesh  used  in  the  PAFEC  analysis  is  presented  in  Figure  4.  The  focus  is  on  the 
critical  region  above  the  crack  in  the  reduced  configuration  shown  in  Figure  2.  Mesh 
refinement  in  this  region  is  necessary  due  to  the  high  stress  gradients  there. 
Optimisation  of  the  angles  used  in  the  refinements  and  the  choice  of  total  mesh  length 
(37mm)  are  discussed  in  appendix  0.  Appendix  0  includes  a  listing  of  the  PAFEC  data 
file  used  in  the  initial  plastic  case.  The  0.02mm  offset  of  the  left  most  nodes  in  the  plate 
relative  to  those  in  the  reinforcement,  was  discussed  earher  with  reference  to  the 
Hjr(0,y)=0  boimdary  condition  across  the  reinforcement. 


(a)  20  14 


18 


DSTOTR-0528 


3.2  Benchmark;  Two-sided  Elastic  Case 

The  FE  analysis  required  an  initial  elastic  calculation  before  each  plastic  run  to 
ascertain  the  load  initiating  plasticity.  The  plastic  run  then  started  just  below  this  and 
increments  took  the  load  into  the  plastic  regime.  Each  elastic  case  was  run  at  the 
reference  load  Fr  (Figure  2):  an  example  is  shown  in  Figure  5  with  the  (elastic)  load 
scaled  up  to  220%  of  Fr,  the  same  as  the  maximum  load  appHed  in  the  corresponding 
plastic  case.  The  stresses  of  particular  interest  are  oii  in  tiie  reinforcement  along  the  top 
surface  and  the  reinforcement-adhesive  interface,  and  oiy  in  the  adhesive,  also  along 
the  interface.  The  stress  components  contributing  to  the  von  Mises  stress,  equations  (2) 
and  (4),  are  also  presented. 


r  .  .  .  .  . . . . _ I _ _ _ _ _ , _ , _ ] 

0  5  10  15  20 


Distance  from  crack,  x  (mm) 


Figure  5.  The  finite  element  stresses  calculated  for  an  elastic  adhesive  at  a  load  of 220%  Fr. 


In  most  cases,  as  demonstrated  in  appendix  C,  the  elastic  cases  gave  results  broadly  in 
tine  with  the  one  dimensional  theory  for  load  transfer  lengths. 


3.3  Two-sided  Plastic  Cases 

Two  cases  were  examined  with  bi-linear  adhesive  stress-strain  responses  (Figure  6). 
The  post-yield  incremental  moduli  were  E-E/10  and  E -E/100  respectively.  These 
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produced  very  similar  results  so  that  "elastic  -  perfectly  plastic"  was  taken  to  be 
represented  by  the  E -E/100  results.  There  was  only  a  small  difference  in  these  two 
cases  when  the  growth  of  the  plastic  zone  length  against  applied  load  (F/Fr)  was 
plotted  in  Figure  9.  The  following  figures  present  the  data  from  calculations  for  the 
two-sided  (symmetric)  joint  case.  The  stresses  of  particular  interest  include  ctxx  in  the 
reinforcement,  both  adjacent  to  the  adhesive  (nodes  8-17)  and  on  the  top  surface  (IS¬ 
IS),  as  well  as  oiy  in  the  adhesive.  The  von  Mises  stress  and  its  decomposition 
according  to  equation  (3)  are  also  plotted  in  Figure  7. 


Figure  6.  Elastic-plastic  stress  behaviour  showing  Young's  modulus  E  and  the  post-yield 
modulus  E'. 

In  this  case,  the  smaUness  of  the  additional  stress  component  oi,  indicates  that  the 
dominant  stress  governing  adhesive  yield  is  axy.  This  explains  the  accuracy  of  the  one¬ 
dimensional  model  (3):  the  adhesive  is  5delding  essentially  one-dimensionally.  The 
axx  stress  along  the  reinforcement  surface  is  interesting  because  it  is  measurable  in  a 
real  specimen  using  strain  gauges  while  the  adhesive  and  reinforcement  stresses  along 
the  adhesive  line  are  not.  If  the  behaviours  at  the  top  and  bottom  of  the  reinforcement 
and  in  the  adhesive  can  be  correlated,  the  onset  or  extent  of  plasticity  may  be 
ascertained  easily  in  practice.  The  reinforcement  surface  stress  peaks  away  from 
directly  above  the  crack  in  the  plate  because  there  is  always  some  bending  of  the 
reinforcement  at  this  point.  This  bending  towards  the  crack  causes  a  lower  (tensile) 
stress  at  the  surface  directly  above  the  crack.  The  other  stresses  in  the  adhesive,  oix  and 
CTyy,  are  relatively  small. 


3.4  One-Sided  Elastic  and  Plastic  Cases 

These  cases  were  tmdertaken  to  investigate  the  magnitude  of  the  peel  stresses  due  to 
bending.  They  used  the  same  mesh  as  the  earlier  calculations,  and  were  derived  by 
replacing  the  %(a:,0)=0  constraint  in  Figure  2  by  Uy  (0.02,0)=0.  They  therefore 
correspond  to  one-sided  repairs  of  plates  half  as  thick  as  for  the  two-sided  repairs.  As 
shown  in  Figure  8,  the  stress  Oyy  in  the  adhesive  is  much  larger  than  in  the  two-sided 
case. 
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t — . . 1 ^ ^ I ^ ^ , , I ^ ^ . 1 

0  5  10  15  20 


Distance  from  crack,  x  (mm) 

Figure  7.  Finite  element  calculated  stresses  for  the  two-sided  plastic  case  with  an  applied  load 
of 220%  Fr  and  post-yield  modulus  E'=E/100. 


t — . — -t— ^ — . — I — ^ . I . . 

0  5  10  15  20 


Distance  from  crack,  x  (mm) 

Figure  8.  Stresses  in  the  one-sided  elastic  case.  The  applied  load  is  again  220%  Fr. 
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The  high  stresses  in  the  adhesive  suggested  that  plastic  yielding  wotild  occur  at  much 
lower  loads  than  for  the  two-sided  case.  This  was  found  as  plastic  yield  would  first 
occur  here  at  only  4%  Fr,  compared  to  the  104%  previously.  To  perform  the  plastic 
calculations,  the  increments  were  reduced  from  10%  load  to  3%  to  ensure  convergence 
at  each  load  increment.  The  load  was  taken  up  to  40%  only.  Clearly,  the  von  Mises 
stress  decomposition  indicates  a  large  contribution  by  non-shear  stresses.  This  case 
needs  further  investigation  and  is  of  interest  because  in  practice  bending  can  occur, 
and  often  only  a  one-sided  repair  is  possible. 

3.5  Discussion  of  Results  to  Date 

The  above  calculations  lead  to  several  useful  derived  quantities,  in  particular  the 
growth  of  plastic  zone  with  increasing  load.  The  one  dimensional  model  leads  to  a 
linear  relationship  as  indicated  earHer  in  Figure  3.  The  growth  in  plastic  zone  length 
with  apphed  load  is  shown  in  Figure  9  for  the  cases  examined.  The  gradient  of  the  line 
fits  that  predicted  by  the  one  dimensional  model  very  well  in  the  two-sided  cases.  This 
reinforces  the  earHer  suggestion  that  the  model  is  accurate  because  the  dominant  stress 
in  the  adhesive  is  oiy.  For  the  one-sided  case,  plastic  yielding  commences  at  a  much 
lower  stress  level  and  spreads  faster  due  to  the  effect  of  the  peel  stress  contributing  to 
the  von  Mises  yield  criterion. 


Figure  9.  Plastic  zone  length  growth  with  applied  load. 
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While  the  plastic  zone  should  not  extend  for  a  large  proportion  of  the  reinforcement 
length,  its  presence  has  the  desirable  effect  of  lowering  the  peak  stress  in  the 
reinforcement.  This  is  because  adhesive  yield  reduces  the  load  being  transferred  into 
the  reinforcement  at  that  point.  Transference  of  this  load  occurs  over  a  longer  section 
of  reinforcement  than  it  would  have  without  plasticity.  The  cr„  stress  distribution 
across  the  reinforcement  is  thus  more  uniform  and  not  so  concentrated  next  to  the 
adhesive.  This  is  particularly  relevant  above  the  crack  in  the  plate.  This  stress  relief 
may  be  shown  (Figure  10)  in  two  ways  as  follows. 

1.  The  first  method  is  by  plotting  the  <Jxx  stress  across  the  reinforcement  at  a:=0  (above 
the  crack). 

2.  Plotting  the  normalised  peak  stress  against  normalised  load  provides  a  more 
graphic  illustration.  The  normalised  peak  stress  is  that  in  the  reinforcement  just 
above  the  crack,  divided  by  the  average  stress  across  the  reinforcement  along  x=0. 
The  normalised  load  is  the  current  load  relative  to  that  needed  to  initiate  plastic 
5delding. 

Figure  10b  shows  a  significant  lowering  of  the  normalised  peak  stress  in  the 
reinforcement,  from  1.6  down  to  1.2  at  220%  Fr.  This  reduction  may  be  an  important 
factor  in  the  choice  of  adhesive  for  this  t3rpe  of  repair.  The  adhesive  thickness  and 
stiffness  will  also  affect  the  stress  concentration.  Calculations  with  eitiier  a  thinner  or 
stiffer  adhesive  layer  produced  higher  peak  stress  values. 


(a) 


(b) 


CU 


o 


c 

o 

e 

o 

o 

c 

’53 

Pi 


Height  above  crack,  y  (mm)  Normalised  load,  F/F 

Figure  10.  Stress  relief  through  jjlastic  yielding,  (a)  The  Uxx  stress  across  the  reinforcement  at 
x=0  for  increasing  loads,  (b)  Normalised  peak  stress  and  surface  stress  against 
normalised  load. 
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4.  Configuration  Variations 

So  far,  this  investigation  has  focussed  on  a  single  repair  configuration,  with  the  various 
thicknesses  and  moduli  chosen  being  liiose  of  representative  repairs.  This  section 
examines  the  effects  of  varying  these  parameters.  In  particular,  the  load  transfer 
lengths  obtained  vary  consistently  with  die  (engineering)  calculated  values  from 
equation  (5)  when  the  various  moduli  are  varied.  The  actual  values  calculated  using; 
(a)  the  formula,  (b)  those  fitted  to  the  exponential  decay  of  the  shear  stress,  and  (c) 
those  derived  from  the  growth  in  plastic  zone  with  applied  load,  are  compared  in 
Figure  11. 


Figure  11.  Variation  in  load  transfer  length,  for  varying  moduli  of  (a)  adhesive,  (b)  plate, 
and  (c)  reinforcement.  In  each  case,  the  modulus  is  expressed  relative  to  that  given 
in  Table  1.  For  example,  E'^reiative=E'^/(1.890GPa). 
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The  load  transfer  lengths  from  the  exponential  decay  and  fitted  to  the  plastic  zone 
growth,  agree  quahtatively  with  the  calailated  values  in  their  behaviour  with 
variations  in  the  moduH.  Quantitative  agreement  among  the  three  values  is  poor.  Table 
2  compares  the  load  transfer  lengths  for  the  initial  case. 


Table  2.  Comparison  between  values  for  the  had  transfer  length  as  determined  by  the  three 
methods,  for  the  initial  two-sided  case. 


Method  of  determination 

Value  (mm) 

Calctilated  by  equation  (5) 

5.881 

FE  result  from  stress  decay 

6.50 

Fitted  to  plastic  zone  growth 

7.82 

The  stress  concentration  above  the  crack,  in  the  elastic  case  prior  to  plastic  yield,  also 
depends  on  the  configuration.  This  is  illustrated  in  Figure  12. 


4.1  Plastic  Zone  Length:  Growth  with  Applied  Load 

As  indicated  earlier,  a  longer  plastic  zone  may  be  beneficial  as  it  leads  to  a  reduced 
stress  concentration  in  the  reinforcement  above  the  crack.  The  zone  should  not  be  so 
long  that  it  extends  for  a  large  proportion  of  the  overlap,  or  else  the  joint  is  hable  to 
failure.  The  one-dimensional  model  predicts  that  tiiis  behaviour  is  hnear  with  applied 
load,  and  for  the  initial  case  shown  in  Figure  13,  this  is  closely  followed  by  the  FE 
results. 

Finding  a  configuration-independent  behaviour  for  this  plastic  zone  growtti  is 
desirable.  To  this  end,  the  one  dimensional  model  suggests  that  plotting  plastic  zone 
length  normahsed  by  the  load  transfer  length,  against  load  normalised  by  the  peld 
load,  should  lead  to  a  "universal  plot".  This  was  indeed  the  case  for  varying  the 
adhesive  yield  stress,  and  was  followed  fairly  well  for  the  other  cases  except  for  the 
stiffer  adhesive  moduli  (Figure  14).  The  problem  with  the  4E'’  case  in  particular,  seems 
to  be  a  delay  in  the  growth  of  the  plastic  zone  with  applied  load  above  the  initial  yield 
load.  It  is  not  clear  why  this  is  so,  but  the  subsequent  growth  rate  agrees  well  with  the 
theory. 
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Varying  adhesive  modulus 


Normalised  load 


Figure  14.  Finite  element  calculated  normalised  growth  of  the  plastic  zone  length  against 
normalised  load  for  the  various  configurational  parameters,  showing  behaviour 
fairly  close  to  theory  in  most  cases. 
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5.  Conclusion 


Stress  reductions  in  the  reinforcement  near  the  crack  region  of  the  order  of  25%  (see  for 
example  Figure  10b)  were  obtained  through  plastic  yielding  of  the  adhesive.  In  the 
two-sided  repair  case,  this  yielding  was  dominated  by  shear  stress,  oiy„  explaining  the 
success  of  the  Hart-Smith  one-dimensional  theory.  The  theoretical  linear  increase  in 
plastic  zone  length  with  applied  load  was  in  good  agreement  with  the  FE  results. 
However  the  observed  load  transfer  length  was  6-18%  longer  than  that  predicted  by 
die  theory. 
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Appendix  A:  Refinements  to  the  Mesh 

The  mesh  for  the  PAFEC  routines  uses  8-point  quadrilateral  elements,  PAFEC  t3q)e 
36210.  These  preferably  have  near-equal  side  lengths  and  angles  near  90°.  The  number 
of  elements  required  is  high  due  to  the  thinness  relative  to  length  of  the  adhesive  layer 
that  has  two  elements  across  the  thickness.  Halving  Ihe  dimensions  of  aU  elements 
thereby  quadrupling  the  total  number  of  elements  and  leading  to  four  across  the 
adhesive,  led  to  results  insignificantly  different  to  those  using  the  mesh  as  presented. 


A.I.  Element  angles 

The  simple  method  of  obtaining  mesh  refinement  near  a  point  of  interest,  jdelds  angles 
of  45°  and  135°  as  shown  in  Figure  15a.  This  can  be  modified  as  shown  in  Figure  15b 
by  constraining  the  points  A  and  B  to  He  on  the  45°  line  and  varying  angles  a,  d  and  e. 
Examining  the  location  of  point  A  first,  its  location  fixes  angle  a.  Angles  h  and  c  are 
then  determined  respectively  as  27Q°-2a  and  45°-a.  Ideally  aU  three  angles  would  be  90° 
but  this  is  not  possible.  The  best  option  is  to  minimise  their  deviations  from  that  by 
minimising 

z  =  2{a-  9oY  +{b-  90Y  +  2{c  -  90)" .  (8) 

The  weightings  arise  because  there  are  two  angles  a  and  c.  The  result  is  c=78.8°, 
b=123.8°  and  c=112.4°. 

A  similar  two-variable  procedure  was  carried  out  with  angles  d  and  e  leading  to 
d=87.4°,  e=71.5°  and/=135°-e=63.5°.  A  preferred  approach  was  to  set  e=f=67.5°  so  that 
the  other  three  angles  at  B  could  be  set  equal  at  75°  and  d=82.5°. 


(a)  (b) 


Figure  15.  (a)  Standard  and  (b)  optimised  meshes  for  refinement  near  a  comer  point.  In  (b),  the 
mesh  has  been  omitted  for  clarity. 
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A.2.  Mesh  length 

The  original  length  of  the  mesh  was  20mm.  With  a  load  transfer  length  arormd  6.5mm 
however,  the  edge  was  affecting  the  results  as  indicated  by  very  different  load  transfer 
lengths  obtained  for  the  elastic  and  elastic-plastic  cases.  Extension  to  37mm  by 
doubling  the  length  (and  number  of  elements)  between  nodes  2  and  15  brought  these 
load  transfer  lengths  into  agreement. 


A.3.  Load  increments 

For  the  plastic  calculations,  PAFEC  begins  with  a  nm  using  an  applied  stress  below 
that  needed  to  cause  plastic  yield.  This  stress  was  calculated  from  the  earlier  elastic 
rtm.  The  load  is  then  incremented  into  plasticity.  The  increments  must  be  small 
enough  to  allow  convergence  and  accuracy  at  each  step.  With  the  von  Mises  yield 
criterion  at  a^=37MFa.,  an  initial  elastic  calculation  indicated  first  plastic  yield  at  62% 
applied  load.  For  plasticity,  load  increments  of  5%,  from  60%  to  90%  produced 
excellent  agreement  with  the  results  using  10%  increments.  The  latter  were  thus  taken 
to  be  sufficiently  accurate  and  increments  calculated  up  to  160%  load.  Using  the  higher 
yield  criterion  of  cry=64.1MPa,  first  peld  occurred  at  104%.  Plastic  calculations 
commenced  from  an  initial  100%  load  up  to  220%,  again  in  10%  increments. 

In  the  one-sided  lap  joint,  10%  increments  led  to  irregularities  in  the  stress  plots.  These 
were  eliminated  using  3%  increments,  but  computing  limitations  restricted  the  number 
of  increments  to  about  12.  Initial  plastic  yield  was  also  much  lower,  imder  10%  load. 
This  case  was  therefore  not  loaded  to  high  levels  (maximum  around  50%)  and  further 
investigation  was  suspended. 
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Appendix  B:  PAFEC  Data  File  for  the  Base  (Plastic) 

Case 


The  following  data  file  listing  is  for  the  base  case  with  plastic  yielding  at  the  von  Mises 
yield  criterion  of  cr'"=64.1MPa.  The  parameters  obtained  in  this  case  are  shown  in  bold 
in  Appendix  C. 


C  DATA  FCR  -  a  PLANE  STRESS  BOOH)  LAP  JOINT  SPECI®! 
C  CPTIMISH)  ICSH  AROID  EM)  WITH  BREAK  IN  PLATE 
C  BASE  CASE,  WU  MISES^A-IIPA 
C 


PLANE.STRAIN 

PLASTIC 

C  C.SLCPE=-(FRACT  OF  E)  TO  APPROACH  STRESS-STRAIN: 

C  0=  INIT  STRAIN  (KXC  HCRIZNTAL),  INF=  INIT  STRESS 
C  («A€  \H?TICALLY  DCMI) 

CC.SL£PE= 

PHASE=1 

PHASE=3 

CPHASE=3 

PHASE=4 


C 

NCDES 

CKI>E.NOeBI  MdS.NOGER  X  Y  Z  THia(NESS 


NCDE.NIICER 

C  PLATE  NCDES 
C  THE  PUTE 

X 

Y 

1 

0.02 

0 

2 

3.12 

0 

3 

0.02 

2.8 

4 

0.27 

2.® 

C  BEUEEN  PLATE  AM)  A)HESIVE 

5 

0.02 

3.1 

6 

0.32 

3.1 

7 

3.12 

3.1 

C  BEILEEN  WHESIVE  AND  PATCH 

8 

0 

3.3 

9 

0.3 

3.3 

10 

3.1 

3.3 

C  PATCH  BCDY  AM)  TCP  SURFACE 

11 

0 

3.6 

12 

0.2 

3.55 

13 

0 

4.57 

14 

3.1 

4.57 

C  POINTS  AT  X=INFINITY 

15 

37.12 

0 

16 

37.12 

3.1 

17 

37.1 

3.3 

18 

37.1 

4.57 

19 

1.27 

3.3 

a  1.27  4.57 

C  VARIABLE  NCDES  TO  CPTIMISE  ANGLES 
CNCDE  X  Y 

C  IN  PLATE  (13,14,1,3,15) 

C  21  0.4281+  0.86S5*T  0 


OCNTRa 

C22 

0.02 

3.1-  1.4142*T 

FUJ..CCNTR0L 

ca 

0.02+ T 

3.1-  T 

TClH?ANCE=iaE-4 

C24 

3.02 

2.7060-  0.8683*7 

C  25  0.02+  1.4142*T  3.1 

C  IN  PATCH  (7,16,8,17,18) 

C  25  1.4142*1  3.3 

CZ7  0  3.3+  1.4142*1 

C  a  T  3.3+  T 

C  29  1.27  3.4672+  0.8685*T 

C  30  0.1672+  0.86£S»T  4.57 

C  SLB-DIVISICNS  OF  BLOCKS  9,11,12 
C  Y(31)=  Y(24),  Y(32)=  Y(29)=  Y(33) 
CWITH  T=  1.8 


C  PHASE=5 

21 

1.99 

0 

PHASE=6 

22 

0.02 

0.55 

PHASE=7 

a 

1.82 

1.3 

C  PHASE=« 

24 

3.12 

1.13 

PHASE=9 

2 

2,56 

3.1 

C  NaN.LI)EAR.T0LERANCE=1  (17#=  DEF) 

C  WITH  T=  0.6 

BASE=4000000 

2 

0.S 

3.3 

CPHASE=10 

27 

0 

4.15 

C  PIGS.STRESS.FILE 

a 

0.6 

3.9 

STOP 

a 

1.27 

3.99 

CCNTRQL.QO 

3) 

0.69 

4.57 

C  BLOCKS  9,11,12  SIB-OIVISICNS 


31 

37.12 

1.13 

32 

3.1 

3.99 

33 

37.1 

3.99 

C  FCR  STRESS  IN  /CIESIVE 

BO 

34 

37.115 

3.150 

2 

37.110 

3.200 

36 

37.106 

3.250 

CBO  NCDES 
C 

PAFBLOCKS 

C  BL0CK.NLraER  TYPE  GRajP.NOBER  ElBBir.TYPE  PROPBITIES 
C  N1  N2  N3  M4  )6  TCPOOGY 
C  REFHS  TO  MATERIAL  AM)  PLATES  AM)  MESH  BELOW 
C  FRCPB?riES:11 
TYPE=1 

ELB©fr.TYPE=36210 

C 

C  NI^-AXIS,  N2  Y-AXIS  TO  DISCRETISE  UOB?  ICSH 
C(REFe®KE)  BL,BR,TL,TR,BC,LO,RC,TC: 

C  B=BOnc«,L=l£FT,R=RIGHT,T=TCP,0=CBnRE 
C  IN  PUTE,  1=LEFT,  feWJACBJT  ADIESIVE,  5=T0  RIGHT 
BLOCK.  GRDLP.  PRCPBTTIES  N1  N2  TOPOLOGY 

1  7  11  1  10  22,23,3,4 

2  7  11  1  3  3,4,5,6 

3  7  11  2  3  4,23,6,2 


2  7  11  1 

3  7  11  2 

C  WIESIVE,  4=LEFT,  5=NEXT 

4  3  12  1 


2  3  4,23,6,2 
XT 

1  4  5,6,8,9 
13  4  6,7,9,10 
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C  PATCH,  frlEFT  WJACEMT,  S^ABCWE 


6 

6 

13 

1 

3 

8,9,11,12 

7 

6 

13 

8 

3 

9,26,12,28 

8 

6 

13 

1 

8 

11,12,27,28 

C  FCR  X=INFINI1Y  TO  RIGHT 

9 

7 

11 

6 

3 

24,31,7,16 

10 

3 

12 

6 

4 

7,16,10,17 

11 

6 

13 

6 

3 

10,17,32,33 

12 

6 

13 

9 

3 

19,10,29,32 

CEXTRA  FCR 

CPTIMISH)  ANGIES 

C  PLATE 

13 

7 

11 

1 

5 

1,21,22,23 

14 

7 

11 

14 

5 

21,2,23,24 

15 

7 

11 

14 

3 

3,24,25,7 

CRATCH 

16 

6 

13 

15 

3 

26,19,28,» 

17 

6 

13 

1 

5 

27,3,13,30 

18 

6 

13 

15 

5 

3,29,30,3 

C  SLB-DIVISICNS  OF  BLOCKS  9 

,12,11=  19,20,21 

19 

7 

11 

6 

5 

2,15,24,31 

20 

6 

13 

9 

5 

29,32,3,14 

21 

6 

13 

6 

5 

32,33,14,18 

CBDPAFBLOCKS 

C 

f€SH 

C  SINGLE  DIGIT  =  N3.  OF  INTERVALS,  LIST=RATIO 
REFBSNCE  SPACING.LIST 
C  X  LEFT  BO  <1,2,4,6,8,  13,17)  NOTE  (BLOCK) 

1  4 

C  X0),  INOUH)  WITH  14  IN  fCSH  13 

2  8,11,15,21,21,21,40,43,44 

C  IN  PLATE  Y(2,3,15)  AM)  PATCH  Y(6,7,13)  INCUDE 
C-6  IN  11,12 

3  4 

C  IN  WfESIVE  Y(4,5,10) 

4  2 

C  IN  PLATE  Y(13,14)  AM)  PATCH  Y(17,18)  INCUDE 
C+3  IN  11,12 

5  2 

C  X  TO  INFINITY  X(9,10,11,19,21) 

6  68 

C  X  TO  INFINITY  X(10),  EXTRA  CVER  ICSH  6  AS  A)TESIVE  THIN 
C  *«*****lJMJSa)  NOW 

7  34 

C  EXTRA  BLOCKS  7,8 
C  X(7)  AM)  Y(8),  INCUDH)  IN  FESH  13 

8  8,11,15,21 

C  X(12)  FILL-IN  PATCH  BLOCK 

9  40,43,44,28,28 

C  REVB1SE  OF  FCSH  2  FCR  X(1) 

10  44,43,40,21,21,21,15,11,8 

C  Y(9)  INCUDES  ICSHES  5,3  ******lJNJSe)  MW,  Y(9)=  TESH 
C  (3) 

11  56,56,49,49,49,49 

C  Y(11,12)  INCUDES  TESTES  3,5  '«*****lJMJSa)  MW,  =TESH 
C  (3) 

12  17,17,17,17,29,29 

C  X(5)  INCUDES  TESTES  (2,14)  AM)  (8,15,9) 

13  8,11,15,21,21,21,40,43,44,28,28 
C  PUTE  X(14,15),  INCUDB)  IN  TESH  13 

14  2 

C  PATCH  X(16,18),  INCUDE)  IN  TESH  13 

15  21,21 
C  EM)  TESH 

C 

TttTERIAL 

C  HATERIAL.NLABB?  E  NU  RO  ALPHA  MJ  K  SH  BULK.H 
C  UNITS  OF  TPa  AM)  inn 


TLATERIAL-TUEER  E  NU 

11  72400  0.33 

C  TD=  (E/2)/  (1-*nj)=  700  TEa  FRCM  PETER  CHALKLEY  FM  73 
CATROCMT 

12  ISO  0.35 

13  207000  0.30 
CBOTKTERIAL 

C 

PLATES.AM).STELLS 

C  PLA1E.MJEai  MATERIAL.TUEER  THIOCTESS  FACING.TWTERIAL 
C  OJTER.LAYER.THICKTESS 

C  RADI  WD2  (LAST  5  COLLTMS  FCR  SAM)WICH  CCNSTRLCTICN) 
PWTE  MATERIAL.TUE®  THICKTESS 
11  11  100 

12  12  100 

13  13  100 

C  END  PLATESJVD.STELLS 
C 

RESTRAINTS 

C  NCDE.NJEER  PLATE  AXIS.MJEER  DIRECTICN  SET 
C  PLANE=0  (NODE  ONLY,  1-3=  X-Z,  4-fr:  LITE  ALCNG  X-Z) 

C  DIRECTICN=  WUl,  1-3=  U(-UZ,  4-6=  ANS(-  ANGZ 
NCDE.TUEB^  PLATE  DIRECTICN 
C  LEFT  EM)  OF  PATOH 
8  1  1 
C  CENTRE  LITE  OF  PLATE 
1  2  2 
C  END  RESTRAINTS 
C 

SLRFACE.FCR.PRESSLRE 

C  LQN).C»SE  PRESSLRE.VAUE  START  FINISH  STE>  NODE  PLATE 

C  AXIS  N1  N2  TS  LIST.OF.NCDES 

PRESSLRE.VAUE  NCDE  PLATE 

C  -100  13  2 

C  STRESS  IN  PLATE:  ALL  3  ARE  E/1000 

-72.4  15  1 

C  STRESS  IN  PATCH 

-207  17  1 

C  STRESS  IN  WTESIVE 

-1.89  34  1 

-1.89  35  1 

-1.89  36  1 

C  EM)  PRESSURE 

C 

INCRBENTAL 

C  G«JSS.POINT.PRINT  DISPLACBOIT .PRINT  NCDAL.STRESS.PRINT 
C  TITE.STEP  SICP.TITE  T«X.STB>.NLFBER  LQW.mSE 
C  SIEP.LIST 

C*.PRINT=  PRINT  FREQ  (DEF=3:  1,4,7,,,), 

C  TITE.STEP,  STCP.T  AND  TWX.STEP.TUEB?  =  CREff  CNLY 
C  PARWETERS,  LQW.CASE  ?,  STEP.LIST=  SIZES  OF  LOAD 
C  INCRBENTS  AS  %  OF  TOTAL  LOAD  CASE(S-  PAIRS  ETC  OF 
C  FRCPCRTICNS  OF  BASIS-SET  LOWS) 

SIH>.LISr 

100,10,10,10,10,10,10,10,10,10,10,10,10 
C  EM)  INCRBENTAL 
C 

PLASriC.T1ATe?IAL 

C  PUSTIC.TWIERIAL.MIEER  TYPE.OF.PLASTIC.FftTffilAL 
C  PRCPER1Y.LIST 

C  PLAS.T1AT.  REFBTRED  BY  YIEUUNG.ELBENTS,  1YPE.OF=  1  FCR 
C  ELASTIC  (IN  TWTERIAL)  AM)  POINr+  POST-YIELD  SLOPES, 

C  2=  STRESS-STRAIN  COORDS, 

C  PRCP.LIST=  CCRRESPCM)ING  DATA  FCR  TYPE.OF 
C  appears  CUT  OF  DATE,  TEH)  DIFFeTBfT 

C  HEWIERS  AS  BELOW  TO  WORK 

PLASTIC.I«TERIAL  YIEU).CRITH?ICN  UNIAXIAL.PRCPeiTIES 

12  1  1 
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C  SO  PLASriC.m'TERIAL 
UmAXIAL.PLflSriC.PRCPB?riES 
UNWXIAL  TYPE  PRCPERTY.LIST 
C  *****note  yield:  shear=37M^,  vm  Mises=64.1 
C  1  1  37,18.9 

1  1  64.1,18.9 

C  BD  lJNIAXIflL.PLASriC.PRCPERTIES 
C 

YIELDING.ELMMTS 

C  PLASTIC.mTERIAL.NLI«a?  OSff.LftJ.NUffiER  START  FINISH 
C  STH>  GROUP.lUeBl  LIST.OF.ELBENTS 
C  PL.WT  REFBS  TO  PL.WT.ND  ABOE  (0=  THAT  OF  EACH 
C  ELBBJT),  CPE3>  FOR  CRE’ 


C  ONLY,  REST=  SPECIFYING  PLASTIC  ElflBfTS  YIELDING 
PLASTIC.MATSIAL.NLMBBI  START  FINISH  STEP 


12 

118 

1 

12 

439 

474 

1 

12 

507 

542 

1 

C  EM)  YIELDING.E1BBITS 
C 

STATE.DeiERMINATICN 

ALQCRITW 

1 

C 

BD.OF.DATA 
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Appendix  C:  Tables  of  Parameters  from  Cases  Studied 


File 

name 

Value 

set 

Calc 

(mm) 

Stress 

(mm) 

Fit 

(mm) 

Calc 

Fy/Fr 

(%) 

Stress 

Fy/Fr 

(%) 

Fit 

Fy/Fr 

(%) 

Fit 

grad 

Varying  adhesive  yield  stress  (Value  in  MPa' 

pll7 

17.32 

5.881 

6.45 

7.17 

26.20 

28.23 

30.5 

0.235 

plong 

37.00 

5.881 

6.40 

8.17 

55.96 

60.3 

67.0 

0.122 

plnew 

64.1 

5.881 

6.50 

Tsz 

96.95 

104.4 

115 

0.068 

pllOO 

100.0 

5.881 

6.40 

7.70 

151.25 

163.0 

179 

0.043 

Varying  adhesive  modulus 

e4 

Ea/4 

11.762 

10.9 

13.18 

193.94 

195.3 

207 

0.0637 

e 

Ea/2 

8.317 

8.95 

9.60 

137.14 

146.5 

150 

0.064 

2e 

2Ea 

4.159 

4.85 

6.14 

68.58 

72.4 

89 

0.069 

4e 

4Ea 

2.941 

3.45 

4.34 

48.49 

49.6 

70 

0.062 

Varying  plate  modulus 

ep2 

Ep/2 

4.740 

5.41 

7.10 

156.32 

162.2 

200 

0.0355 

2ep 

2Ep 

6.883 

7.43 

8.26 

56.75 

60.6 

63.5 

0.130 

Varying  reinforcement  modulus 

er2 

Er/2 

4.867 

5.50 

6.57 

80.25 

87.65 

98 

0.067 

2er 

2Er 

6.704 

7.45 

8.58 

110.54 

116.6 

128 

0.067 

Varying  adhesive  thickness 

t 

2tA 

8.317 

8.85 

10.12 

137.14 

138.85  151 

0.067 

C.l.  Notes  to,  and  definitions  of  terms  in,  the  above  table. 


The  row  in  bold  is  the  "base  case":  all  parameters  were  varied  from  this  case  as 
indicated. 

Calc:  calculated  using  configuration  and  one  dimensional  model  formulae, 

stress:  obtained  from  oiy  stress  data  (Figure  5, 7). 

Fit:  linear  fit  to  plastic  zone  growth  plot  (Figure  9). 

(mm):  Load  transfer  length  in  mm 

Fy/pR  (%):  5deld  load  as  a  percentage  of  the  reference  load  (far  field  E/1000  stresses 
applied:  Figmre  2). 

Fit  grad:  gradient  of  the  linear  fit  to  the  plot  of  plastic  zone  length  in  mm,  against 
applied  load. 
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